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I* ilSIocgcTiow 



Histcrically , advanced datsction and surprise attacks 
have teen an area cf interest to the military community. 
Knowledge about the effects of atmospheric boundary layer 
(A3L) inversions on elect rcmagnetic signals is well known 
and dccumented. Advance knowledge of inversion changes can 
have a direct effect on military planning conducted ty all 
services. 

Kith technological advances, our weaponry has become 
highly sophisticated (e.g., the development cf "smart" 
bombs) and dependant on electromagnetic propagation for 
guidance ccntrol. The guidance components for these weapons 
depend cn many different wavelengths, ranging from the 
visible to the microwave, dependent on the specific use of 
the weapon. The weapons depend on the transmission and 
reception cf these signals to locate and follow targets, 
Cifficulties arise when atmospheric conditions are such that 
signal performance is distorted by refraction. Only those 
areas which affect propagation of radar signals will be 
discussed in this thesis. Other propagation types, such as 
optical, and the aticspheric parameters affecting th<=m will 
not be addressed. 

Signal propagation is dependent on the index cf refrac- 
tion which is a function of temperature and humidity. 
Eefracticn cf a signal occurs as that signal moves through 
vertical gradients cf temperature and humidity. Atmospheric 
conditions which enhance refraction are inversions, which 
causs strong gradients of temperature and humidity, or rapid 
changes cf humidizy near the surface. Occurrence, height 
and intensity of inversions are strongly dependent on local 
and synoptic atmospheric conditions. The importance of 




.1 




predictir.g propagaticr is obvicus but to do sc with accuracy 
first laquirss accurate orsdiction of the mer ecrclogical 




prsdicticn of changes in 



■Of obvious i rpcrtr. r.c i~ oh? 
the temperature and humidity 



profiles, both of which are highly variable and can 
completaly change their vertical profile within a matter of 
hours. Other parameters such as vertical velocity and 
surface wind speeds will be discussed later. 

The ability to forecast inversions and how they may vary 
directly influences the effectiveness of modern weaponry. 
It is readily apparent that improvements in this forecasting 
ability is cf paramount importance in making military plans 
which involve choosing the type weapon best suited for a 
particular m-issicn. As long as there are functions depen- 
dent cn signal transmissions there will remain a need to 
increase the ability to forecast the atmospheric parameters 
which ccrtribute to propagation. 

The model used to investigate inversion heights is a 
marine atmospheric boundary layer (HABL) model which has 
been modified for use at the Naval Postgraduate School 
(NPS) . The factors incorporated into this model include 
potential temperature and specific humidity within the 
well-mixed region and their gradients at the base of the 
inversion. Use of these parameters make this model appli- 
cable to studying inversion changes. Actual discussion of 
the model will be provided in Chapter II. 

Another area where advance knowledge of ABL inversion 
changes would be useful is in single-station forecasting. 
The effects of chances in the inversion are found not only 
to effect electronic transmissions but also to indicate 
possible changes to currenx weather conditions. Therefore, 
any icprcvements in cur ability to predict these changes 
will be important to land-based and sea-based operations, 
since imprcvement in this area cf prediction also will be 
reflected in the daily weather forecast. 
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This th=sis complst^s an analysis on four tima paricds 
cf data obtained frcni Vandenberg 'aFB, California. These 
data periods each cover a different season and vary frctr tvic 
to eight days in length. In conducting this exaainatior. 
vertical velocities fcr Vandenberg AFB were calculated from 
their tavinsonde data. Each days data were run through the 
DPS HAEL prediction model to establish the forecast trend 
for that period. Jesuits from the model's forecasts were 
compared to data frci later reports for correlation and to 
determine accuracy cf the model. Successive runs of the 
model were made until each data set has been completed and 
determination as to how well the model predicts inversion 
heights in each 24 hcur forecast. 

One cf the time periods of data chosen for study is, by 
choice, coincident with a data period used by Gleason 
(1982). His study was directed at effectiveness of the 
model for more of a pure marina location and fcr this reason 
data from San Nicolas Island was chosen. This study uses 
data from Vandenberg AFB since it is a coastal station and 
will have somewhat different surface conditions due to the 
continental influence. See Fig. 1.1 for a geographical 
orientation of the area of interest. Once the Vandenberg 
AFB data had been studied, the San Nicolas Island results 
were used to make comparisons. The results cf this compar- 
ison are of interest since they will indicate whether cr not 
there are potential problem with the ilABL model when used in 
a regime strongly influenced by land as compared to a regime 
closely resembling marine ccndirions. 
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II. HOD EL DESC RIPT ION 



1. EICKGRCDUD 

The MAEL model in use at the Naval Postgraduate School 
(NPS) is a zero-order, two layer, integrated mixed layer 
model (Davidson et ai. 1984). This is a modified version of 
the model developed by Stage and Businger (1981). The two 
layers consist of the well-mixed, turbulent boundary layer 
and the relatively ncn-turbulent free atmosphere. 
Separating these layers is the transition zone known as the 
inversion. The term, ”ze ro- order" , applies to this model 
because the inversion is considered a jump from one value to 
another instead of having a finioe gradient. Thus, we 
assume that, at the inversion, the eguivalent potential 
temperature ^nd the total specific humidity (g ) 

undergo zero-thickness jumps. 

Total specific humidity and equivalent potential temper- 
ature are conserved quantities in pseudo-adiabatic processes 
and, for this reason, are assumed well-mixed within the 
boundary layer. They are given by: 

'it = '^1 (2.1) 



where 



is the water vapor mixing 



liquid water mixing ratio; and 



ratio and q is the 



0 

e 



+ L / (C q ) 

V p V 



( 2 . 2 ) 



where 9 is the potential temperature, 
for vaporization of water, and C is 
constant pressure of air. a third 



is the latent heat 
the specific heat at 
fundamental dependent 
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Figure 2. 1 Typical potential temperature and 
moisture profiles. 

variable is the mean mixed layer depth (h) . Typical 
profiles of these parameters are shown in Fig. 2.1. As can 

he seer in Fig. 2. 1, both 0^ and g^ are assumed zc have 
constant values throughout the mixed layer. Above this is 
the interfacial region, or inversion, where the jump occurs. 

One cf the more appealing features about this model is 
that it runs from input parameters gleaned from rourine 
metscrclcgical and cceanographic dana. These parameters 
are: surface layer humidity, temperature and winds, 

vertical gradients cf temperature and humidity above xhe 
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itixsd lay'=r, and s=a-surfacs temperature (SST) . Readily 
available rawinscnde cfaservaticns can be used zo detarjiine 
ail input parameters excepr SST. SST values can be cbrair.ed 
from satellite int erpretati cn cr direc:; measurements. 

The predictive equations for the three dependent vari- 
ables are given by: 

dq /dt = h“^( (Wq') + w Aq^) 

t. u e t 

d9^/dt = h (W9') + R - R + w A6 ) (2.3) 

c u c G © 

dh/dt — w + w 
D e 

where d/dt represents the time rate of change of the vari- 
able following a column of air. 

surface turbulent fluxes of 0^ and g^, while the A terms are 
the magnitude of the jumps at the inversion. Other quanti- 
ties introduced here are: 1) w - the entrainment rate; 

2) Wj^ - the large-scale vertical velocity at the inversion 
height (initialized using the subsidence rates described 
later in this chapter); 3) Eg - the rate of radiative heat 
loss per unit area near the cloud top; and 4) E^ - the 
rate of radiative heat gain per unit area at the cloud base. 
(Haggerty, 1983) 

B. NES MODEL PARTICDIABS 

The use of a well-mixed layer, along with other davelcp- 
ments, made modifications to the model described by Stage 
and Eusinger (1981) necessary. The NPS integrated HABL 
differs from others in the specifics of surface flux, radia- 
tive flux and the zero gradient assumption. The following 
discussions are brief since the references (Davidson et al. 
1984) are quire complete. 
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St iiicsp heric structu r<= 



lb's profile ci a well-oiixed atmospheric variable, -j) , 
is r e c r e £ 31! tel as 



4> (z) 
({> (z) 



<j) + z 

ms 9m 



mh 



r . (z-h) 

9 



Z<h 

Z>h 



(2.U) 



where m refers to the mixed layer, s to the surface value, h 
is the height of the mixed layer, Z is rhe height, is the 
jump in across the inversicn zone, is the gradient 
above the inversion, and is the gradient in the mixed 
layer (assumed tc be zero) . 



2 . Dyn ami c s within the mixed layer 

The model predicts the time evolution of the mixed 
layer values of g^, 0^, and h (see Egn. 2.3). At initiali- 
zation, the lapse rates -for these guantiries are held 
constant while larger, synoptic scale observations provide 
the input values of the boundary layer wind for the entire 
forecast period. 

Surface fluxes of sensible heat and water vapcr are 
determined through an application of the bulk aerodynamic 
formula (Davidson et al. 1978). These fluxes are driven by 
the bulk air/sea differences and are expressed as 






e 0 



m 






(2.5) 



where Cg is the ten meter, stability dependent, drag coeffi- 
cient, u* is the friction velocity, and 9^ and g^ are the 
respective sea-surface values. 
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Par ametsriza-ticn for the ent rairnien-t rate is thar 
given fcy Stage and Evjsinger (1981) * This paratrie t= rizat ion 
allows closure of the system of equations and determines the 
time evolution of the inversion height. It also brings the 

turbulent kinetic energy (TKE) budget into the formulation 
since entrainment of heat flux across the stable inversion 
xequires TKE, that is, work is done by this process. Here 
the assumption used for closure states that the dissipation 
rate of TKE is a fixed fraction (1-A) of the production 
rate, where A is the entrainment coefficient, set equal to 
0.2. Within the TKE equation are terms which provide for 
production and dissipation, therefore it is assumed that 
both occur within the mixed layer, rather than only where 
there is negative buoyancy flux. These are associated 
through the use of the turbulent scaling velocity for the 
mixed layer, w* , and the mixed layer depth (h) for the 

scaling length. 

3 ♦ Bad iati v a flux cons ider at ions 

When considering radiation flux and its effects on 
the mixed layer there are several points which must be 
understood. Distinction must be made between long-wave 

fluxes for clear or cloudy cases, and short-wave fluxes. 

a. Long-wave radiative flux 

Dnder cloudless sky conditions, the net long- 
wave radiative flux is determined using the water vapor and 
temperature profiles. The flux is calculated by inte- 
grating, from the surface (Z=0) to the top of the mixed 
layer (Z=h) , the flux emissivity profile where the emis- 

sivity is determined from the Stef an-Boltzmann law (Fleagle 
and Eusinger 1980) . The aforementioned modification has 
also been applied. 
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The presence of clcuds necessitates the calcula- 
tion of the net radiative flux at both the cloud top ana th'-=> 
cloud base. These calcul aticns are also made using the 
Stef an-Ecltzmann law. When the cloud top flux is being 
determined, xhe input temperatures are the cloud rep temper- 
ature and the effective radiarive sky temperature. Fcr the 
cloud base the input temperatures are rhe rem pe natures at 
the cloud base and at the sea surface. It is also important 
to note that, in this case, the flux divergence between the 
sea surface and the cloud base is neglected. 

b. Short-wave radiative flux 

The MABL modal calculations of short-wave radia- 
tion take into consideration both direct and diffuse compo- 
nents. Eoth scattering and absorption contribute to 
short-wave extinction. The scattering contribution tc the 
short-wave radiative flux is determined through use of the 
delta-Eddington method, where the mixed layer scatrering is 
due rc atmospheric particles, cloud droplets in the cloudy 
oase and saa-salt aerosols in rhe clear sky case, 
flbsorpticn is primarily the result of water vapor and is 
taken into account when the incident flux at the top of the 
mixed layer is determined. 

C. MCEEI AFPLICATIOS 

All model runs for this thesis ware done such that rhe 
situation approximated to a real forecasr run. That is, 
the model was initially set tc run with a single sounding 
and an estimate of -C.0005 m/s for the subsidence. Further 
runs were made as each new sounding could be added. This 
presents a realistic representation of what one could expect 
from the model for the Vandenberg AFB area. 
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Ereparation for running the model requires that several 
points or items bs considered. These consist of identifying 
th-a synoptic pattern and initialisation parameters; '^'ili be 
discussed below. 

1 . Ida nti f y S yn op tic P attern 

Eefcre checking on input data for the model the 
overall synoptic pattern must be checked. To do this cne 
utilizes all available meteorological data (e.g, surface and 
upper air synoptic analyses, as well as satellite, rawin- 
sonde and individual surface observations) . The desired 
synoptic situation fcr optimum model predictability is to 
have the area of interest dominated by a high pressure 
syste IT . 

2* Ide nti fy Ini tializa tion P arameters 

The next step in preparing to make the model runs is 
to deteriine the input parameters. Potential temperature 
and moisture values fcr the mixed layer, height and strength 
ci the jump, surface wind, sea-surface temperature, vertical 
velocity and the lapse rate above the jump are the necessary 
parameters and, of these, a good determination of the jump 
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represent the true sounding or serious error can result. 

The vertical velocity or, in our case, the subsi- 
dence rate for the forecast starx time must be determined. 
There are several means for calculating vertical velocity 
but none, tc date, provide the accuracy desired. With xhis 
in mind, the method which has proved to be most accurate was 
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used, Ihis method is the Q-method which is an expression 
for vertical velccitj derived from the vertical integration 
cf the ircisture budget aquation and is given by: 

A (q + l) (dh/dt) - h(d(q + l)/dt) + (Wq’) 

w = ii (2.6) 

A(q+1) + 3(h/2) 

where g and 1 are vapor and liquid water contents, and 6, 
their vertical gradient {Gleason, 1982). The computer 
program used to calculate the vertical velocity uses current 
data as well as data within the previous 24 hours. The use 
cf previcus rawinsor.de data is allowed since rhe Q-method 
equaticn includes total change with respect to time (d /dt) 
for tha inversion height and the total moisture. The 
vertical velocities output frcm the program are for each 
data period input. These values are rhen averaged to get 




figure 2.2 Sample of digitized rawinsonde. 
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the rat€ fcr ths er.tirs period with the final sutsicer.ce 
rate calculated teinc the one used in the model. The only 
difficulty encountered with this method was when the subsi- 
dence inversion height was increasing at a rate sufiicisnn 
to yield a positive average. When this occurred the model 

input rate was set tc zero. Table I gives the subsidence 

rates (ic/s) used for the model runs. 



TABLE I 

Sutsidence Rates for Hodel Runs 



CATE 

July 



Dece irber 



April 





GUT 


SOESIDNECE 


DATE 


GMT 


SUBSIDENCE 


25 


00 


-0.C05 


August 29 


00 


-0.005 




1 2 


-0.003948 




1 2 


-0. 0C08393 


26 


00 


-0.003359 


30 


00 


0.00 




12 


-0.000826 5 




12 


-0. 0002389 


27 


00 


-0.01489 


31 


00 


-0.0004 126 




1 2 


-0.001553 




12 


0.0 


28 


00 


-0.002856 


September 1 


00 


0.0 








2 


00 


-0.005934 










12 


-0.001676 


13 


00 


-0.C05 


3 


00 


-0.000082 




1 2 


-0.005813 




12 


0.0 


14 


00 


-0.002856 


4 


00 


-0.0006649 




12 


0.0 




1 2 


-0.0002713 








5 


00 


-0.0005856 










12 


-0.001524 


12 


00 


-0.C05 


6 


00 


0.0 




1 2 


-0.002019 








14 


00 


0.0 










1 2 


-0.000365 1 








15 


00 


O.C 










12 


0.0 








16 


00 


-0.C02208 










1 2 


-0.002988 
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III. SJROFTIC P&TTEBS DISCUSSION 



fcur data sets were used for evaluation of model 
validity. These sets include a summer period (25 through 27 
July 1977), a winter period (13 and 14 December 1977), a 
spring period (13 through 17 April 1981), and a fall period 
(29 Acg through 5 September 1981). Large scale analyses 
were for as close to 1200 Greenwich Mean Tims (GMT) as 
possible. In most cases the satellite pictures will he at 
1745 GMT since these pictures best display the area of 
interest . 

This discussion will be of the overall weather picture 
for the eastern North Pacific ocean and the western United 
States. Included will be conditions over the area as a 

whole and unique or differing conditions at Vandenberg AFB 
itself. Satellite pictures, surface weather analyses, and 
upper air analyses (SCO millibars) will be used to determine 
the synoptic situation of each period. 

i. SOHMIE CASE, PEHIOD I (25-27 JULY 1977) 

The surface and 500 millibar (mb) synoptic patterns for 
this period appear in Figs. 3.2 through 3.4. During this 
period the coastal region was dominated by a high pressure 
system over the North Pacific Ocean with the thermal 
troughing extending from Mexico northward into the high 
desert region of Arircna and southeastern California. Such 
a synoptic pattern usually produces weak onshore flow along 
the southern Califorria coastal regions. The upper air 
analyses (500 mb) showed no strong westerly flew over 
California during the entire period. 
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Durirg this pericd ~hs Pacific high pressure system 
moved next heas* ward and on 27 July# weakened as a 
stationary front developed eff the northwest coast. Also of 
importance cn 27 July was ths development of a low pressure 
center west of Los Angeles, California Even with this 
development the preceminate flew remained onshore in the 
Vanderberg AFB region. 

The rawinsonde data revealed a subsidence type inversion 
throughout the pericd, see Fig. 3.1. Thera was alsc a 
diurnal variation in the inversion base. Each day the base 
ranged from the surface at 1200 GMT to a height above 650 
feet at COOC GMT. An example of this diurnal change also 
appears in Fig. 3. 1 . 

The satellite pictures for this period appear in Figs. 
3.5 through 3.7. Orly patchy areas of lew clouds along the 
Califerria coast appear on the satellite pictures. These 
areas were in the San Francisco bay area and the coastline 
south cf Lcs Angeles. The only ether area cf low clcuds was 
two tc three degrees longitude west of Vandenberg AFB. 
There was nc significanr changes to the coastal locations of 
the lew clcuds while the area over the water had some dissi- 
pation and was even further offshore. The Vandenberg AFB 
vicinity remained under clear skies during the entire 
period. 
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(a) 0000 GMT Sounding 




POT. TEMP. (C) 0 10 20 30 40 



(t) 1200 GHT Sounding 

Figure 3.1 Bxanple of Sunnier Case rawinsonde and 
diurnal change in inversion. 
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Piqnre 3,2 Surfaca <a) and 500 mb (b) analyses 
1200 G3T 25 July 1977. 
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(a) 




Piqure 3.3 Surface (a) and 500 mb (b) analyse 
1200 GHT 26 July 1977. 
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(a) 




Fiqtixe 3.4 Surface (a) and 500 nb (b) analyses 
12C0 GBT 27 July 1977. 



29 




Figure 3.5 Satellite picture for 25 July 1977, 
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ligure 3.6 



Satellite 



picture for 26 



July 1977. 
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figure 3„ 7 Satellite picture for 27 July 1977. 
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SISIIH CAS2, PEB2CD II |13-14 DEC2HBEE 1977) 



Ih€ ger.sral surface and 500 mb analyses appear in 
figures 3.9 end 3.10 A difference between this syncopic 
patten and that fcund in the other periods is in this 
patrern the desired onshore flow exists but it dees not 
arise from the strong Pacific high and thermal rroughing tut 
from a high pressure cell that had separated frem the main 
Eacific cell and meved in over central California Ihe 
strongest westerly flow aloft was found to the north of 
Califcrrla with the Polar jet meving across the Washington - 
CregcE region. Westerly flow over Vandenberg AFB was 
stronger than in other periods. 

Ihe primary aspect is that on 13 December a high pres- 
sure center moved over central California with the large 
high pressure system still maintained over the North Pacific 
Cesan. Weak thermal troughing was also extending northward 
into California The overall result of these surface 
patterns ^as weak onshore flew in the Vandenberg AFB and San 
Nicolas Island regions. The only other significant feature 
in the 13 Eecemfcer analysis was the frontal system moving 
into the Pacific Northwest. 

Cn 14 December the high pressure center dissipated 
leaving the Pacific high pressure system and ridging across 
central California. Weak thermal troughing continued in the 
southern California area. Offshore flow also continued to 
dominate the southern California coastline. The frontal 
system from the 13th moved eastward and a trailing warm 
front moved to northern California. 

Fig. 3.8 is an example of the vertical soundings for 
this period. The rawinsonde data showed subsidence type 
inversiers on both days. Most important is that on 13 
December at 0000 GMI the inversion was very weak and not 
well structured. Ihe inversion strengthened and remained 
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aloft until 1200 GMT cn 14 December when it was observed as 
a surface based inversion. At this rime the inversion 
reached the maximum strength for rhe period. 

The satellite pictures for this period appear in 7igs- 
3.11 and 3.12. Clear skies in the Vandenberg AFB area with 
a frontal system moving intc the northern California regions 
appear in the satellite picture for 13 December. There is 
very little indicaticn of low clouds along the southern 
California coastline. By 14 December the front is well into 
northern California with increased cloudiness across central 
Califernia. There has also been an increase in the amount 
of low clouds along the southern California coast. Patchy 
stratus follows the coast south of Los Angeles and a large 
area now sits off the Vandenberg AFB coast. 




POT. TEMP. CO 0 10 20 30 -40 



Figure 3.8 Example of Sinter Case rawinsonde. 
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(a) 




(b) 

Piqure 3.9 Scrface (a) and 500 mb (b) analyses 
1200 GHT 13 December 1977. 
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Piqure 3.10 Surface (a) and 500 ab (b) analyses 
1200 GHT 14 Deceaber 1977. 
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Figure 3.11 Satellite picture for 13 December 1977 
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Figure 3.12 



Satellite 



picture 



for 14 December 1977. 
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SERIBG CASE, PERIOD III (13-16 APRIL 1981) 



Ihe synoptic patterns fcr this period appear in Figs. 
3.14 through 3.17. The dominating feature was a high pres- 
sure system extending nort heasrward from rhe Norxh Pacific 
Ocean regions across Washington into Canada. Thernal 
troughinc dominated the high desert regions of California 
and extended into northern California. The prsdorninate 
southern California coastal flow, generated by these 

patterns, was onshore. 500 mb analyses shew light and vari- 
able winds with tendencies toward formation of a low pres- 
sure center through the period. The strongest westerly flow 
remained well north of the Vandenberg AFB area. 

Cn 13 April the Eacific high was strong but showing a 
tendency to weaken as the longwave pattern extended into 
Canada. As the longwave high pattern shifted to the central 
United States rhe Pacific high weakened. The thermal trough 
in scuthern California continued to dominate but also had 
some weakening. 

Ey 15 April the Pacific high system had begun to 
strengthen with some ridging across central California. The 
high continued intensificat ion through the remainder of the 
period. Cn 16 April the ridging extended into the 
Washington area with the Pacific high pressure center begin- 
ning to icve northeastward. Throughout this portion of the 
period the thermal trough showed very little variation. The 
only ether key occurrence from the 15th to rhe 16th was the 
movement of a cold front across the Pacific Northwest inxo 
Nevada. 

Fig. 3.13 is an example of the rawinsonde soundings for 
this period. Rawinsonde reports during this period had 
subsidence inversions with much daily variation. At the 
beginning of the period the inversion was weak and aloft. 
Eiurnal variations intensified the inversion but also 
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brought it to the surface. This pattern conxiaued through 
the 14th but on the 15th the inversion intensified at JCOO 
GMT and remained aloft until 1200 GMT on the 16th. During 
the latter half of the period there was some fluctuation in 
the intensity but over all the inversion remained strong. 

These satellite pictures appear in Figs. 3.18 through 
3-21. On 13 April the satellite showed Vandenberg AFE under 
clear skies while the coastal region south of there was 
dominated by low clouds. ty 14 April the area of coverage 
by low clouds increased and Vandenberg AF3 was within that 
area. low clouds continued to dominate the Vandenberg AFB 
area with seme dissipation beginning on 15 April. 16 Apr 
showed crly thin coverage of low clouds over the Vandenberg 
AFB area with the denser coverage south of the area. 




Figure 3,13 Example of rawinsonde for Spring Case. 
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(a) 



(b) 

Figure 3.14 Surface la) aad 500 mb (b) analyses 
1200 GHT 13 Jpril 1977. 
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(a) 




Pigure 3.15 Scrfac^ (a) and 500 nb (b) analyses 
12C0 GMT 14 April 1977. 
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(b) 

liqtire 3.16 Surface la) and 500 mb (b) analyses 
1200 GST 15 April 1977. 
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(a) 




(b) 

Fiqtirs 3.17 Surface (a) and 500 mb (b) analyses 
1200 GHT 16 april 1977. 
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ligcre 3.18 Satellite picture for 13 April 1981 
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figure 3. 19 



Satellite picture for 14 April 1981. 
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figure 3,20 



Satellite picture for 15 April 1981. 
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Figure 3.21 



Satellite 



picture 



for 16 April 



1981 . 
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B. PILL CASE, PERIOD IV (29 ADGOST - 5 SEPTEMBER 1981) 

■Ih^ surface (a) and 50 0 mt (b) analyses for Psricd 17 
appear ir Figs. 3.23 through 3.30. The predominate surface 
feature during this sxtire period was rhe presence of a high 
pressure system in the Eastern portions of the North Pacific 
Ccsan. A thermal low and trough system dominated the 
southern California rsgion. The effect of these two surface 
features was to produce weak onshore flow in the area of 
concern. Conclusions drawn from the 500 mb analyses show 
that, except for the last day of the period, the maximum 
westerly flow remained north of California. Only on 5 

September did the Polar jet extend toward the regions of 
interest . 

The surface analysis for 29 August (Fig 3.23) shows the 
western Pacific Ocean dominated by a high pressure system 
with a lew pressure center located off the Washington coast. 
A thermal low extended northwesxward from Mexico into 
Calif erria with troughing extending into northern 
Califernia. The southern end of a cold front had also moved 
into the northern California area. 

Ey 30 August the cold front had moved into southern 
Califerria and the thermal low had changed its orientation. 
The Pacific high pressure system showed weak intensification 
with movement toward the northeast and had pushed inland by 
the 31st. The 31st also shewed that the cold front was 
beginning to leave California but with another frontal 
system moving out cf the Gulf of Alaska. Figures 3.24 
through 3.25 depict these occurrences. 

The 1 September analysis revealed a weakening of the 
high pressure and a cold front approaching Califernia. A 
thermal trough was well established over California hut very 
diffuse on 2 September as the cold front moved through 
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southern California. The high pressure system again inren- 
sifisd and showed nc marked weakening throughcun the 
remainder cf the period. The only significant event after 2 
Septenifcar was the entrance of a rapidly moving cc_i front 
into California. This scenario is shown in Figs 3.26 
through 3.30. 

A saiple of the rawinsondes for this period appears in 
Fig. 2.22. The rawinscndes for this period have a very good 
subsidence type inversion throughout the period. In each 
rawinscnde release ths inversion was found ahcve 1000 feet 
with decreasing relative humidity through the inversion 
zone. 

The extent and development of the low clouds appear in 
the pictures in Figs. 3.31 and 3.38. The satellite pictures 
revealed a buildup cf coastal low clouds (stratus and/or 
fog) during the period. Beginning on 29 August the clouds 
tended to he thin with some dissipation during the day. 
Further into the period the stratus thickens and by ths end 
of the pericd there is a dense area of low clouds off the 
central and southern California coastal regions. 
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SPEC. HUM. (gykg) 5 10 15 20 25 

POT. TEMP. (C) 0 10 20 30 40 



Figxire 3. 22 



Bawinsonde example for Fall Case 





Jiquxe 3.23 Serf ace (a) and 500 «b (b) analyses 

12C0 GHT 29 august 1977. 
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(a) 




(b) 






and 500 ab 
august 1977. 



(b) 



analyses 
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ligure 3.24 



Surface 
1200 GHT 




(a) 




liqtre 3.25 Scrface (a) and 500 mb (b) analyses 
1200 GMT 31 august 1977. 
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(a) 




(b) 



Tigure 3.26 



Sorface (a) aiid 500 mb (b) 
1200 GHT 1 Septeaber 1977. 



analyses 
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liqure 3.27 Surface (a) and 500 mb (b) analyses 
1200 GHT 2 September 1977. 
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(b) 



Figure 3.28 



Surface (a) and 500 mb (b) 
1200 GUT 3 September 1977. 



analyses 
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(a) 




(b) 



rigure 3.29 



Surface (a) and 500 mb (b) analyses 
1200 GHT 4 September 1977. 
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(a) 




(b) 

ligure 3.30 Surface (a) and 500 mb (b) analyses 
1200 GHT 5 September 1977. 
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Figure 3,31 



Satellite picture for 29 August 1981. 
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Figni® 3.32 



Sat€llits picture for 30 August 1981 
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Pigare 3. 33 



Satellite 



picture for 31 



August . 
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Figure 3.34 Satellite picture for 1 September 1981. 
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Figar€ 3.35 Satellite picture for 2 September 1981. 
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Pigar€ 3.36 



Satellite pictnre for 3 September 1981. 



65 









Figure 3.37 



Satellite picture for 4 September 1981. 
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Figure 3.38 



Satellite picture for 5 September 1981. 
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IV. RESOITS 



flcd^l forecasts 'were made for every 30 minutes for 24 
hours following the initialization time. The forecasts 
presented in the following tables show only the 12 and 24 
hour values since they correspond to the verif ica ricn tines. 
The structure of the tables are such that the first fore- 
cast, indicated by cclumn heading FCST1, was based on the 
first observation for initialization and wirhin that column 
are the forecasts generared. For example, if the first 
observation was for 1/0000 GMT then the forecasts from -hat 

observation would be for 1/1200 GMT and 2/0000 GMT. These 

forecas-cs would be in column FCST1 . The second forecast, 
column heading FCST2, was initialized with the second obser- 
vation (e.g., 1/1200 GMT) and within that column are the 

forecasts generated (2/0000 GMT and 2/1200 GMT). Each table 
continues in this manner until all forecasts for a period 
have been given. 

Boot mean square (RMS) scores are given following each 
forecast except the last. No EMS for the last forecast were 
given as there were no observations for comparison. Also 

note that the 12 hour RMS value simply represents the 

difference between observed and forecasted values. 

The values given in the tables are for the parameters: 
inversion height, lifting condensation level, mixed layer 
potential temperature and mixed layer specific humidiry. 
Units for these are meters for the inversion height and the 
lifting condensation level, degrees Centigrade for potential 
temperature, and g/kc for speoific humidity. Also note that 
subsidences, used in discussing these results, are given in 
Table I. 
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m£CCS3I0N OF EOCEl RESDLTS FOR PERIOD I (SOMMER) 






This period «as dcminated by clear skies. The forecast 
lifting condensation level values (see Table III) only fore- 
casted clcuds during FCST2. The forecasted inversion 

heights. Table II, reveal that the model was either very 
close or quite in error relative to the observed height. 
Tables IV and V reveal that, throughout the period, there 
were alsc problems associated with potenrial temperature and 



TABLE II 

July Inversion Results 



DATE 


GKI 


CBS 


FCST1 


RMS 


FCST2 


RMS 


FCST3 


RUS 




25 


00 


356 


- 


- 


- 


- 


- 


- 






12 


157 


195 


38.0 


- 


- 


- 


- 




26 


00 


277 


1 08 


173.2 


92 


185.0 


- 


- 






12 


199 


- 


- 


112 


204.4 


165 


34.0 




21 


00 


206 


— 


— 


— 


— 


100 


111.3 




DATE 


GMT 


CBS' 


FCSTa 


RMS 


FCST5 


RMS 


FCST6 


RilS 


FCST7 


21 


CO 


206 


172 


34. 0 


— 


— 


— 


- 






12 


120 


146 


42. 8 


154 


34.0 


- 


- 


- 


28 


CC 


199 


- 


- 


1 16 


89.7 


71 


128.0 


- 




12 


- 


- 


- 


- 


- 


42 


- 


1 10 


29 


CO 


- 


- 




- 


- 


- 


- 


62 



TABLE III 

July Lifting Condensation Levels 



DATE 


GMT 


OBS 


FCST1 


FCST2 


FCST3 


FCST4 


25 


00 


633 


- 


- 


- 


— 




12 


146 


299 


- 


— 


- 


26 


00 


111 9 


199 


42 


- 


— 




12 


217 


- 


0 


519 


- 


27 


00 


206 


- 


- 


309 


185 




1 2 


452 


— 


— 


— 


168 


DATE 


GMT 


OBS 


FCST5 


FCST6 


FC3T7 




27 


1 2 


452 


605 








28 


00 


126 1 


410 


1 84 


- 






12 


- 


- 


176 


586 





00 - - - 366 



6 9 



29 



TABLE 17 







July 


Eotential 


Temperatures 




DATE 


GMT 


OBS 


FCST1 


ECST2 


FCST3 


FCSTU 


25 


00 


16.2 


- 


- 


- 


• 




12 


11. 6 


17.4 


- 


- 


- 


26 


00 


17 . 4 


18.2 


11.6 


- 


- 




12 


12. 6 


- 


8.5 


17.8 


— 


21 


00 


20. 9 


- 


- 


18.4 


15, 2 




12 


14. 1 


— 


— 


— 


16 . 7 


DATE 


GMT 


055 


FCST5 


ECST6 


FCST7 




27 


12 


14. 1 


19.8 


- 


- 




28 


00 

12 


19.5 


19.6 


17. 1 
18.3 


18.9 




29 


00 


- 


- 


- 


19.3 











TABLE V 










July Moisture Values 




DATE 


GMT 


OBS 


FCST1 


FCST2 


FCST3 


FCST4 


25 


00 


8. 2 


- 


- 


- 


- 




12 


7.9 


6.3 


- 


- 


- 


26 


00 


6 . 8 


5.4 


8.3 


- 


- 




1 2 


8. 1 


- 


7.2 


9.7 


- 


27 


00 


9.3 


- 


- 


11.3 


9,8 




12 


7.9 


- 


— 


— 


10-8 


DATE 


GMT 


OBS 


FCST5 


FCST6 


FCST7 


j 


27 


1 2 


7.9 


10.6 


- 


- 




28 


00 


7.3 


1 1.6 


11.1 


- 






12 


- 


- 


12.0 


10.1 




29 


00 


- 


- 


- 


11.6 





specific humidity forecasts. These problems will be covered 
during the following discussions. 

Examination of the inversion height results reveals seme 
good and seme poor height forecasts at the 12 hour point but 
mostly bad forecasts at 24 hours. Those forecasts made with 
average subsidence rates (-0.003 to -0.005 ra/s) had fair to 
good results at 12 hours but had poor results at 2U hours. 
The best height forecasts (Table II column FCST4) were with 
the weakest subsidence rate while the second best (the fifth 
forecast) was with the strongest subsidence rats. Also of 
interest is that forecasts with the strongest subsidence 
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also hac tha lifting con den 
inversicn. Ancthar iuiportant 
except fcr the forecast with 
best forecasts were for the 1 
1200 GHI (a morning period fcr 
Fctertial temperature, s 
(with and without clcuds) and 
in explaining what was cccurr 
sion heights for this peric 
described separately as no its 
sion forecast. 



ation level well aocve cn= 
feature in Table II is that, 
the weakest subsidsr.c*, the 
hour periods initialized at 
the California region) . 
ecific humidity, radiaticn 
entrainment will be examined 
ng with the predicted inver- 
. Each parameter will be 
possiole impact on the invsr- 



1) potential temperature - Under clear skies one 
would normally expect nightly cooling over land. 
The observed values shown in Table IV reveal that 
this occurred yet the forecasts do not reflect this 
fact. Nocturnal ceding strengthens the jump and 
when coupled with subsidence lowers the inversion as 
seen in Fig. 4.1. All forecasts were lowering the 
inversion due to good subsidence rates for initiali- 
zation. Daytime heating acts against the subsi- 
dence, weakens the jump, and results in lifting of 
the inversion. This diurnal event was not reflected 
in the forecasts. 

2) specific humidity - After the first two fore- 

casts, specific humidity foracasr values were only 
fair. There were no forecast values which were good 
since the errors ranged from just under 1 g/kg up to 
3.3 g/kg. The primary trend in the forecasts was to 
increase specific humidity. This prediction of too 
large a specific humidity caused the model to 
predict even lower lifting condensation levels 
which, in turn, enhances the poor inversion 

pre diet ions. 



71 



3) radiaticE - This factor must be considered nncer 
conditions of clear as well as cloudy skies. 

a) cloudy ~ In the second forecast the rodel 
predicted clouds to fora and extend to the 
surface. The clouds radiate long wave energy to 
space so the potential temperatures were fcrecasr 
too low. Inversion predictions, while decreasing 
then increasing during ohis forecast, did not 
reflect the diurnal change predominate with clear 
skies . 

t) clear - Radiation computations here could 
have been incorrectly influenced by the large 
specific humidity forecasts. This accompanied 
with good subsidence rates caused all forecasts 
to lower the inversion. 

4) entrainment rates - these influence the time 
rate of chance of specific humidity and potential 
temperature as well as inversion height. In several 
cf the forecasts the entrainment rats was zero so it 
was internally set to a positive, yet minimal, 
value. This could have affected the other fore- 
casted values but it is believed the effect was less 
than others. 

fi definite diurnal pattern can be seen (see Fig. 4. 1) in 
the cbserved inversicn heights. Generally there were 
increased heights in the transition from morning to evening 
and decreases fcr the opposite Transition. This pattern 
coincides with the diurnal pattern of tha potential tempera- 
tures which reflect radiative daytime heating and nightly 
cooling. Surface winds also displayed the influence of 
daily heating with evening reports having the strongest wind 
speeds and with steady northwesterly directions (an onshore 
wind). Figure 4.1 also depicts the 12 and 24 hour forecasts 
and it is easily seen than the diurnal pattern shown by the 
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height (meters) 




Figa^e 4,1 G 
Observed an 




73 



cbserv€d data is 



net shc'rfn by the forecasts. Only the 
Eornirg tc evening 12 hour forecast indicate a cerrset 
trend. However close inspection of the forecasts indicate 
that at these times the model inputs for subsidence were 
strong thereby influencing the forecast. 

These effects all suggest problems with the radiation 
calculations threugheut this summer period. A reasen for 
this error is the fact that Vandenberg AFB is a shoreline 
station and the model is primarily designed for over water 
use. Iherefora, land heating and cooling fay solar radia- 
tion, even through a stratus sky, require adjustments to 
apply the model at a shoreline station. Vertical velocity 
and entrainment rates also affected forecasted MABL proper- 
ties through influences on the rate of change of key parame- 
ters as well as the inversion heights (see Sgn. 2.3). 



B. BISCCSSION OF MODEL EESDLTS FOR PERIOD II (WINTER) 

Eesults in Table VI generally suggest that during this 
season the model fails to forecast variations in the height 
of the inversion. Tables VII through IX compare forecast 
and observed lifting condensation levels, potential tempera- 
tures and specific humidity values. 



TABLE ?I 

December Inversion Results 



DATE 


GMT 


CBS 


FCST 1 


1 3 


CO 


164 


- 




12 


171 


45 


14 


00 


284 


14 




12 


297 


- 


15 


00 


— 


- 




12 


- 


- 



RMS FCST2 RMS 
126.0 

210.7 40 244.0 

11 265.8 



FCST3 RMS FCST4 

187 110.0 

124 



7 4 



492 

623 



I 



I 

I 

I 



TABLE VII 

Dec’eabar lifting Condensation Levels 



^ s 


GiiT 


OSS 


FCST1 


fCST2 


FCST3 


FCS ! 4 


13“ 


0 0 


1083 


- 


- 


- 


- 




1 2 


2055 


123 


- 


- 


- 


14 


00 


1056 


22 


936 


- 


- 




12 


1 74 


- 


338 


511 


- 


15 


00 


— 


- 


- 


255 


0 




12 


- 


- 


- 


- 


0 



TABLE VIII 

Deceaber Potential Temperatures 



EiTE 


GMT 


OBS 


FCST1 


FCST2 


FCST3 


FCST4 


13 


00 


17. 2 


- 


- 


- 


- 




1 2 


15. 8 


14.6 


- 


- 


- 


14 


00 


15.7 


14. 1 


15.4 


- 


- 




12 


13.0 


- 


15.2 


14.9 


- 


15 


00 


- 


— 


— 


14.8 


8.2 




12 


- 


- 


- 


- 


5.4 



TABLE IX 

December Hoisture Values 



CATE 


GMT 


OBS 


FCST1 


FCST2 


FCST3 


13 


00 


8.9 


- 


- 


- 




12 


3. 6 


9.7 


- 


- 


14 


00 


6,3 


9.9 


6.6 


- 




12 


8, 3 


- 


9.0 


8.0 


15 


00 


- 


- 


- 


9.2 




12 


- 


- 


-- 


- 



FCST4 

7,3 

7.0 



Vandenberg AFB remained under clear skies in this 
yet the forecasts were for rapid formation of 
Forecast specific humidity values improved in the 
forecast but were still not acceptable for the 
period. The potential temperature forecasts were 
the 12 hour point. 



period 
clouds, 
second 
e nt ir e 
best at 
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Ihe prirnary causes for the poor parf or manes c£ the modal 
in this case appear tc bs vertical vslocicy and sntraintrant. 
Ihe radiation and cloud caicalacions centribansd cc the 
error, tut to a lesser extent. Entrainment was a ksy factor 
during the first three forecasts because null values were 
obtained in computation, sc ninimal positive values were 
assigred. The radiation segment calculations are strongly 
influenced by the presence cf clcuds and on the last day the 
model forecast clouds. With zero vertical velocity, the 
model forecast rapid lifting of the inversion. With no 
ebssrvatiens it is hard to evaluate those forecasts but the 
trend was for increasing heights although not as rapidly as 
predicted. 

Forecast and observed heights for this period appear in 
Figure 4.2. The diurnal variations noted in the summer case 
are net present in this case. It is believed that the large 
errors in the inversion height forecasts are due to two 
factors. First, the vertical velocities with 'which the runs 
were initialized were large which lead to forecasts of 
lowering heights in all but the final run. In reality, the 
heights shewed a continual increase. The vertical velocity 
used in the last run was such that rapid lifting of the 
inversion was forecast. Second, the surface winds which 
stayed strong even through the morning hours allowed mere 
mixing in the Icwer layers thereby increasing the depth of 
the mixed layer. This is especially noticable in the 
13/1200 GMT surface wind (360 degrees at 2.1 m/s) . 

At Vandenberg AFE a small chain cf mountains exists 
due-north of the base. with this northerly flow, there 
would be an increase in the local turbulsnce in the lower 
levels. Therefore, the poor height forecasts in this case 
is associated with the land influence which is net taken 
into acccunt by the nrcdel. 
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C. DISCCSSION OF SODEI EESOLTS FOR PERIOD III (SPBIRG) 

Tatl€s X through XIII list th? fcracasred values for the 
spring case. This period was cloud free the firs- day and 
had lew clouds the second. Ihe forecasts did not include 
cloud formation until the third day. Even rhough the cloud 
forecasts were in error it does not appear that the cloud 
calculations had any real effect on the poor inversion 
height predictions. 

The vertical velocities and entrainment rates were 



TABLE X 

April Inversion Results 



EAIE GMT 


OBS 


FCST1 


RMS 


FCST2 RMS 


FCST3 


RMS 


13 00 


14 1 


- 


- 


- - 


- 


- 


12 


136 


31 


10 5.0 


- - 


- 


- 


14 00 


28S 


9 


299.0 


75 214.0 


- 


- 


12 


150 


- 


- 


40 240.6 


293 


143.0 


15 00 


469 


- 


— 


— - 


298 


222.9 


DATE GMT 


OBS 


FCST4 


R MS 


FCST5 RMS 


FCST6 


RMS 


15 00 


469 


140 


32 S.O 


- - 


- 


- 


12 


509 


128 


503.4 


624 115.0 


- 


- 


16 00 


42 1 


- 


- 


835 429.7 


6 84 


263.0 


12 


291 


- 


— 


— — 


867 


663.2 


DATE GMT 


OBS 


FCST7 


RMS 


FCST8 






16 12 


291 


340 


4 9.0 


— 






17 00 


- 


440 


- 


323 






12 








335 






different than 


those for other 


period. 


Entrainment rates 


were 2 erc fer 


the 


first 


and third forecast 


s, and the first 


three heurs 


of the fi 


fth forecast and 


we re 


assigned a 


aininal positive value of .001. 


Hence, the 


inve 


rsion height 


was subsidence 


conticll 


ed for 


these times 


• 


The moisture 


forecasts were 


good since values agree to 


within 2 g/kg of 


the ctserved value 


a nd i 


n s ome 


cases even 


less 


than 1 g/kg. 


Each of these 


forecasts 


we re 


initia lized 


with 


the evening 



sounding in which lifting condensation levels indicated no 
clouds. 
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TaBLi XI 

April Lifting Condensation Levels 



LJTE 


GMT 


0 B3 


FCST1 


FCST2 


FCST3 


FCST4 


13 


00 


705 


- 


- 


- 


- 




12 


3 70 


103 


- 


- 


- 


14 


00 


989 


59 


274 


- 


- 




12 


305 


- 


191 


919 


- 


15 


00 


506 


- 


- 


861 


234 




12 


445 


— 


— 


— 


1S0 


E3TE 


GMT 


OBS 


FCST5 


FCST& 


FCST7 


FCST8 


15 


12 


445 


315 


- 


- 


— 


16 


00 


623 


U23 


398 


- 


- 




12 


256 


- 


303 


297 


- 


17 


00 


- 


- 


- 


206 


78 




12 


- 


- 


- 


- 


0 









TABLE 


XII 










April 


Potential Teaperatures 




BATE 


GMT 


OBS 


FCST1 


FCST2 


FCST3 


FCST4 


13 


00 


13. 6 


- 


- 


- 


- 




12 


7.3 


12.8 


- 


- 


- 


14 


00 


15.7 


12.9 


10.7 


- 


- 




12 


8.2 


- 


12.5 


15.1 


- 


15 


00 


13. 2 


- 


- 


14.7 


10.3 




1 2 


10.6 


— 


— 


— 


11.4 


BATE 


GMT. 


OBS 


FCST5 


FCST6 


FCST7 


FCST8 


15 


1 2 


10.6 


9.9 


- 


- 


- 


16 


00 


12. 4 


9.6 


8.7 


- 


- 




12 


8.3 


- 


7. 1 


12.3 


- 


17 


00 


- 


- 


- 


9.8 


5.6 




12 


- 


- 


- 


- 


2.7 



Close examination reveals that the potential temperature 
forecasts did not reflect the observed diurnal variaticn. 
The 12-hcur forecast generally showed cooling or heatir.a as 
it should but not sufficiently to agree with observed 
results. The 24-hour forecast tended -o reverse the diurnal 
trend and to have coding instead of heating cr vice-versa. 

It is believed that the land again influenced the 
observed inversion heights. Figure 4.3 shows the cbserved 
and forecasted inversion height values for this case. Close 
examination reveals that there appears to be two problems 
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TABLE ZIII 

Afxil Hoisture Values 



DATZ 


G.:iT 


OBS 


FCST1 




FCST2 


FCST3 


FCST4 




13 


00 


6.6 


- 




- 


- 


- 






12 


5.2 


8.9 




- 


- 


- 




14 


00 


6.5 


9.0 




6.9 


— 


- 






12 


5. 7 


- 




8.2 


6.5 


- 




15 


00 


7.2 


— 




— 


6.6 


6 . 8 






12 


6. 2 


— 




— 


— 


7,6 




DATE 


GMT 


OBS 


FCST5 




FCST6 


PCST7 


FCST8 




15 


12 


6.2 


6.4 




— 


- 






16 


00 


6.4 


5.9 




5. 6 


- 


- 






12 


5.9 


— 




5.3 


7.6 


- 




17 


00 


- 


- 




- 


6.7 


5.4 






12 




■ 






■ 


5,2 




which are 


a combine 


tion of 




the problems dis 


cussed for 


the 


summer and 


winter ca 


sss . 


Ih 


3 £ r r s 't. 


part of 


the period 


was 


strongly influenced 


by diurnal fluctuations 


and the fo 


re- 


casted values were sinilar 


tc 


those in the summer case. 


The 


winter case 


was 


inf 1 


uenced 


by 


su bsid 


ence and 


surface winds. 


In this case, it is 


belisv ed 


the nor 


therly winds caused 


the 


increase in 


the 


observed inversion 


at 15/0000 GMT due 


to 



tcpcgraphic induced large mechanical mixing. The vertical 
velocity allowed the icodel to decrease the heights for that 
time. 



E. EISCESSION OF HODil RESULTS FOR PERIOD IV (FALL) 

Tables XIV - XVII list the observed and forecasted 
values. As before predicted inversion heights for this case 
appear tc be no better than those for previous cases. The 
formation cf clouds was forecasted accurately since ‘ the 
nodal predicted cloud formation at a time when low clouds 
occurred in the area. However, forecasted lifting condensa- 
tion levels were not in agreement with those observed. 
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Figure 4.3 Graphical Depiction pf Spring Case 
Observea and Porecastea laversxon Hsignts. 

8 1 



TABLE XI? 

August-Septembcr Inversion Results 



EATf 




OBS 


FCST1 


R MS 


FCST2 


RMS 


FCST3 


P s 


29 


CO 


287 


- 


- 


- 


- 


- 


- 




12 


297 


137 


16 0. 0 


- 


- 


- 


- 


30 


00 


365 


67 


336.2 


524 


159.0 


- 


- 




12 


3 1 1 


- 


- 


6 38 


363.6 


383 


72.0 


3 1 


00 


391 


- 


— 


— 


— 


4 96 


127,3 


CATE 


GKT 


OBS 


t?c'^T4 
‘ 512 


RMS 


FCST5 


RMS 


FCST6 


RilS 


3 1 


00 


39 1 


12 1.0 


- 


- 


- 


— 




12 


553 


666 


16 5.6 


574 


21 .0 


- 


- 


1 


00 


322 


- 


- 


755 


433.5 


819 


497,0 




12 


M 


• 


— 


— 


— 


1057 


— 


CATE 


G«T 


OBS 


FCST7 


R MS 


FCST8 


RMS 


FCST9 


RMS 


1 


12 


a 


363 


- 


— 


— 


- 


— 


2 


00 


356 


545 


18S.0 


- 


- 


- 


- 




12 


356 


- 


- 


273 


83.0 


- 


- 


3 


00 


4 1 1 


- 


- 


225 


244.2 


449 


38,0 




12 


469 


- 


— 


— 




5 09 


55,2 


CATE 


GMT 


OBS 


FCST10 


RMS 


FCST11 


RMS 


FCST1 2 


RMS 


3 


12 


469 


491 


2 2.0 


- 


- 


- 


- 


4 


00 


422 


714 


292.8 


684 


262.0 


- 


• 




12 


529 


- 


- 


8 54 


417,5 


4 54 


75.0 


5 


00 


424 


— 


— 


— 


— 


562 


157,1 


CATE 


GMT 


OBS 


FCST13 


R MS 


FCST 14 


RMS 


FCST15 




5 


00 


4 24 


■ 675 


2510 


- 


- 


- 






12 


429 


823 


467. 


2 404 


25 


- 




6 


00 


- 


- 


- 


496 


- 


4 98 






12 




- 


- 


- 


- 


561 





Ihe forecasted specific humidity values have an inter- 
esting trend (shewn in Table XVII) . Forecasts initialized 
from evening data «ere good but those initialized from 
morning data were bad. When initialized with evening data, 
the predicted specific humidity values were representative 
cf the ctserved values. Obviously, pcor forecasts of 
specific humidity or any other key parameter lead to pcor 
prediction of inversion heights. 

Ectential temperatures, in all but two runs, were fore- 
cast to always decrease (cod) with time. The two morning 
to evening forecast cases did not reflect sufficient daytime 
heating, even with clouds present. To evaluate the pcor 
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TABLE XV 



Augast-S apta 


aber Lift! 


ng Condensation Levels 


CATE 


G.^T 


OBS 


FCST1 


FCST2 


FCST3 


FC S 1 4 


29 


00 


73 1 


- 


- 


- 


- 




12 


1 96 


503 


- 


- 


- 


iC 


00 


520 


234 


227 


- 


- 




12 


206 


- 


0 


78 


— 


31 


00 


406 


- 


- 


0 


297 




12 


364 


— 




— 


206 


CATE 


GMT 


OBS 


FCST5 


FCST6 


FCST7 


FCST8 


3 1 


12 


36 4 


105 


- 


• 


— 


1 


00 


467 


227 


613 


- 


- 




12 


M 


- 


580 


101 


- 


2 


00 


331 


- 


- 


150 


— 




12 


205 


~ 


- 


- 


0 


3 


00 


577 


— 


— 


— 


0 


DATE 


GMT 


OBS 


FCST9 


ECST10 


FCST1 1 


FCST12 


3 


00 


577 


237 


- 


- 


- 




12 


2 84 


122 


325 


- 


- 


H 


00 


496 


- 


452 


424 


— 




1 2 


452 


~ 


- 


341 


0 


5 


00 


525 


— 


— 


— 


0 


DATE 


GMT 


OBS 


FCST1 3 


FCST14 


FCST15 




5 


00 


525 


374 


- 


- 






12 


310 


263 


394 


- 




6 


00 


- 


- 


186 


194 






12 








0 




per f ormancs 


in potential temperature 


forecasrs, the subsi- 


dence rates 


, specific humidit 


y and d 


iurnal 


radiaticn varia- 


tions were 


com pa 


r 0 cl • 


When the 


forecast was 


for a morning to 


evening fer 


iod , 


excluding the 


two pre 


viously 


mentioned fore- 


casts, ccn 


tinued coding was 


predicted when heating should 


have been p 


redic ted . 


This result was indep 


endent of subsi- 


dence rates 


or presence of clouds. 


Therefore, for morning 


to evenirg 


forecasts 


, radiati 


on effe 


czs wer 


e not sufficient 


to introduce the 


needed dayti 


me heating. 


From evening to 


morning the major 


problem was too 


much 


cooling. This 


happened in 


five of 


the forecasts with ano 


ther five fore- 


casts shewing good 


ceding rates but 


bad forecasts for the 


24-hctr poi 


nt since 


the previ 


ous 12- 


hour fo 


recast was bad. 


The geed ceding 


rates occurred when 


the s 


ubsidence was at 


its largest 


compared to the 


excessive lowering forecasts. 


This wculd 


in dicat 


€ possible probl 


9ms in 


the radiation 
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TABLE XVI 



dugus 


t-ssptember Totential 


Temperat 


ures 


E ATE 


GMT 


OBS 


FCST1 


FCST2 


FCST3 


FCSTU 


29 


00 


17. £ 


- 


- 


- 


- 




12 


13.5 


16.6 


- 


- 


- 


30 


00 


16. 6 


15.5 


10.4 


- 


- 




12 


13.8 


- 


7.6 


13. 7 


- 


31 


00 


15. £ 


- 


- 


1K5 


9.8 




12 


14.2 


— 


— 


— 


7.6 


EATE 


GMT 


OBS 


FCST5 


FCST6 


FCST7 


FCST8 


21 


12 


14.2 


12.3 


— 


- 


— 


1 


00 


17. 6 


13.2 


12.1 


- 


- 




12 


M 


- 


10.3 


13, 8 


- 


o 


00 


16. 1 


- 


- 


12,5 


- 




1 2 


12.8 


— 


- 




8 . 6 


3 


00 


15. 4 


— 


— 


— 


9.7 


DATE 


GMT 


OBS 


FCST9 


fCSTIO 


FCST1 1 


FCST12 


3 


00 


15.4 


9.9 


— 


- 


— 




12 


13.1 


7.8 


12.2 


— 


— 


4 


00 


16. C 


- 


10.8 


10. 1 


- 




12 


14.1 


- 


- 


7.9 


12.2 


5 


00 


16. £ 


— 


— 


— 


11.2 


DATE 


GMT 


OBS 


FCST1 3 


FCST14 


FCST15 




5 


00 


16. € 


11.6 


- 


- 






12 


13.8 


9.7 


16.1 


- 




6 


00 


— 


- 


12.8 


11.2 






12 








8.7 




which 


would 


be cool 


ing the 


pczenti 


al tsmp 



for this ti^e of day. 

Jigure 4.4 shows the observed and forecasted inversions 
for this period. Fey factors affecting inversion heights 
during this period were diurnal heating variations, surface 
winds, and clouds. The temperatures showed a definite 
diurnal pattern throughout the period as did the lifting 
condensation level. The inversion displayed a diurnal trend 
only part of the time, 29/1200 GMT through 31/0000 GMT and 
3/0000 GMT through the remainder of the period. 

The low subsidence values contributed greatly to the 
high inversion height values forecasted during this period. 
Of the remaining forecasts, those with heights forecast to 
decrease were evening to morning forecasts with sutsicence 
occurring. The observed heights increased (except for the 

2/OOOC GMT value which remained constant) due to strong 
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TABLE XVII 





A' 


ugust- 


September 


Moista 


re 7alues 


DATE 


GMT 


OBS 


PCST1 


FCST2 


FCST3 


FCSI4 


2S 


00 


8.5 


- 


- 


- 


- 




12 


3.7 


9.0 


- 


- 


- 


3C 


00 


9 . 1 


9.7 


6.9 


- 


- 




12 


8. 8 


- 


6.6 


9.3 


- 


31 


00 


9.0 


- 


- 


9.3 


6.4 




12 


8. 3 


— 


— 


— 


5.8 


DATE 


GMT 


OBS 


FCST5 


fCST6 


FCST7 


FCST8 


:• 1 


12 


8.3 


8,4 


— 


— 


— 


1 


00 


9.9 


8.3 


6 • 3 


- 


- 




12 


M 


- 


5.7 


9.3 


- 


2 


00 


9.6 


- 


- 


8.3 


- 




12 


8.2 


- 


- 


- 


9 . 1 


3 


00 


8.0 


— 


— 


— 


8.9 


CATE 


GMT 


OBS 


FCST9 


FCST10 


FCST1 1 


FCST12 


3 


00 


8.0 


6.7 


- 


- 


- 




12 


8. 1 


6.1 


7.4 


- 


- 


4 


00 


8.7 


- 


6.3 


6.1 


- 




12 


7.9 


- 


- 


5.5 


9.0 


5 


00 


9.0 


— 


— 


— 


9 . 1 


CATE 


GMT 


OBS 


FCST1 3 


FCST14 


FCST15 




5 


00 


9.0 


6.9 


— 


— 






12 


8. 3 


6.4 


9.3 


- 




6 


00 


- 


— 


8,3 


7.5 






12 


- 


- 


- 


7.0 





winds (wsaker and more westerly at 2/0000 GMT) which lifted 
the inversicn due to lixing . 

At 31/1200 GMT an observed height decrease occurred for 
the morning to evening forecast time but the forecasts 
showed rises. To explain this requires consideration of all 
three factors mentioned above. The observed height values 
decreased due to the increase in cloud coverage and 
continued onshore flow whereas the normal situation is 
offshore flow due to diurnal changes. The radiation compu- 
tations also appear to have been responsible for errors in 
the forecasts for this period. Radiarion was not the only 
problem since erroneous specific humidity and vertical 
velocities contributed. There was land influence on the 
parameters involved in height predicrions, as in the other 
cases . 
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Figure 4.4 
Observed 



Graphical Depiction of 
and Forecasted Inversion 



Fall Case 
Heights. 
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E. SOMMIBl AND GENEEAI COaBENlS 

It is easily seen by the iargs root mean squara sccras 
3 d h 9 ‘ii ? d 9 T9 nc^ s bs'dw^^r. dor9css*tz Ti d cbstrirv^d 

trends, depicted graphically, thar the model did ncr yislc a 
good prediction of inversion heights at Vandsnberg AF3. The 
model also did not yield very good forecasts of the ether 
parameters. 

The discussions present an explanation of the prcbatls 
components where errors occurred. There are certain prime 
factors which causa the greatest problem. The largest error 
producer is the fact that Vanbenberg AFB is a land base 
which strongly experiences diurnal variations in temperature 
and winds. The primary predictive values of the model 
(subsidence, entrainment, radiation and presence of clcuds) 
are also influenced by the land effect. Daytime heating 
which increases temperature and, hence, dissipates clcuds is 
not accernted for in the MAEL model. 

Topography and certain wind directions increase the 
turbulence in the mixed layer. With mountains north as well 
as south of Vandenberc AFB, this condition can result when 
the offshore flow does not set up during the night and 
strong northerly or southerly flow continues. This leads to 
deepening of the mixed layer and can overcome the effects of 
subsidence. The result of this would be to lift instead of 
lower the inversion. 

It is concluded that this study yielded conclusive 
evidence that the MAEl model is not appropriate for use at a 
coastal station because, even though ccasral, land influ- 
ences play a significant role. The results from the study 
by Gleasen (1982) reveals definite contrasts with this one. 
He had geed results for San Nicolas Island and was able to 
show that the Q-methcd of calculating vertical velocities 
was best. 
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soaa aRY ahd concldsions 

Accurate predicticns of variations in the height of 
atmospheric inversions are a necessity in proper usage of 
modern weaponry. This is because of the impacr changes in 

inversion heights have on the ducting of weapons guidance 
sysrem or of tracking sysrem signals. Atmospheric parame- 
ters are continually changing and those parameters which 
govern inversions can easily change their entire profile in 
a matter of hours. Forecasts cf changes in inversion prop- 
erties trust start vith forecasts of those parameters, 
namely, temperature and specific humidity. To do this the 
marine atmospheric boundary layer model used at the Naval 
Eostgraduate School was utilized. 

Data periods which covered each season were chosen for 
the Vandenberg AFB area. All four cases were dominated by a 
high pressure system over the eastern North Pacific Ocean. 
Clear as well as cloudy days were considered to allow full 
range testing for a station strongly influenced by land 
conditions. Data for these periods were from rawinsonde 
observations, satellite information, and surface and 500 mb 
synoptic analyses. Model input parameters were from the 
rawinsonde observations and satellite derived sea-surface 
temperatures. 

Eesults from each model run were compared with observed 
values tc determine accuracies of the model. It was found 
that the model did not yield accurate predictions during 
most cf the cases. Final interpretations of these inaccura- 
cies always point to strong land influence. Diurnal varia- 
tions influenced by the land are much greater than 
variations over water and are not reflected in the model 
calculations. 
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Ih€ €nd result cf this study is that this MAE! nodal 
dees net '<«ork for Vancenberg AFB. It is not currently valid 
for stations strongly influenced by land induced diurnal 
variations. The redevelopment cf the model to treat such 
factors is required to properly handle the land influences. 
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